Tail Probabilities and Partial Moments for Quadratic Forms in Multivariate Generalized Hyperbolic Random Vectors

Countless test statistics can be written as quadratic forms in certain random vectors, or ratios thereof. Consequently, their distribution has received considerable attention in the literature. Except for a few special cases, no closed-form expression for the cdf exists, and one resorts to numerical methods. Traditionally the problem is analyzed under the assumption of joint Gaussianity; the algorithm that is usually employed is that of Imhof (1961). The present manuscript generalizes this result to the case of multivariate generalized hyperbolic (MGHyp) random vectors. The MGHyp is a very flexible distribution which nests, among others, the multivariate t, Laplace, and variance gamma distributions. An expression for the first partial moment is also obtained, which plays a vital role in financial risk management. The proof involves a generalization of the classic inversion formula due to Gil-Pelaez (1951). Two applications are considered: first, the finite-sample distribution of the 2SLS estimator of a structural parameter. Second, the Value at Risk and Expected Shortfall of a quadratic portfolio with heavy-tailed risk factors.